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1. Introduction 

A persistent discrepancy of about 3.2 a between the measured [1] and the theoretically calcu- 
lated [2] anomalous magnetic moment of the muon, a^, has motivated several lattice calculations 
of the leading order hadronic contribution, Aa^ [3-7]. This contribution is the dominant source 
of error in the theoretical estimate of and presents both a challenge and opportunity for lattice 
QCD calculations. Here we extend our previous work [6] to address most of the systematic errors 
in an effort to establish what is required for a precise calculation of Aa^. A detailed analysis will 
be presented elsewhere, so these proceedings will focus on the contributions of the vector mesons. 

2. Leading order hadronic contribution to 

The leading order hadronic contribution to the anomalous magnetic moment of the muon can 
be calculated from the vacuum polarization tensor using the following expression from [3]. 

f°° 1 

Aa M = a 2 / dQ 2 — w (Q 2 / m j)U R (Q 2 ) (2.1) 
Jo Q z 

The renormalized vacuum polarization function is TIr(Q 2 ) = H(Q 2 ) — n(0) and the weight func- 
tion is given in [3]. The relevant lattice details for this calculation can be found in [6]. 

The hadronic correction to an is normally estimated using a dispersion relation and experi- 
mental measurements of the cross section o(e + e~ — > hadrons). A recent review [2] provides a 
list of results ranging from Aa^ = 684.6 (6.4) • 1(T 1() to 711.0 (5.8) • 1(T 10 . The variation in the 
experimental estimates is beyond our current precision, so in the following we will use just the 
estimate of 691.0 (5.3) • 10~ 10 from [8]. More important than the total value for Aa^ are the indi- 
vidual contributions. The low energy region is dominated by the vector mesons: p, ft) and 0. The 
contributions from each are given in [8] as 

AaJ =501.1 (3.0) -10~ 10 Aa M | ft) = 37.0 (1.1) - 10 10 Aa^| =34.42 (0.9) -10~ 10 . 

Taken together, these three contributions already account for about 83% of the total Aa^. Veri- 
fying the contributions from the vector mesons is clearly a good first step towards calculating the 
complete Aa^. 

3. Vector mesons 

As discussed in previous lattice calculations of Aa^, the low Q 2 region of Hr(Q 2 ) requires 
an interpolation and extrapolation in order to perform the integration in Eq. |2.1| . Currently this 
requires introducing some model assumptions. Given the significant contributions of the lightest 
vector mesons, it is reasonable to incorporate their contributions into the Q 2 shape of the model. 
First we establish what is known regarding the electromagnetic decays of the p, ft) and (p. Then 
we discuss the additional model assumptions needed to calculate the contributions of the vector 
mesons to Aa u . 
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3.1 Vector meson decay constants 

The electromagnetic coupling of a vector meson is defined by 

(n|/£ m (0)|V,/?,£) = mygv, em £n(p) fv,em = m v gV,em ■ 



The value of gv,em can be determined through the decay V — > e + e . A straightforward exercise 
gives the partial width for this decay as 

An 

y 



r(V -> e + e ) = — a 2 m gl em m v (\+2m 2 Jml)J l-Am 2 Jm 2 . 



Using the latest values from the PDG [9], we find 

g p , em = 0.20174 (86) g(0 , em = 0.05863 (98) g^ em = 0.07473 (120) . 
It is common to represent these decay constants in an isospin basis defined as 

j'= Q = (uy ll u + dy ll d)/\ f 2 j'= l = (uy^u -dy^d) / '\fl J^sy^s. 
Assuming m u = m^, is a pure ss state and ft) has no ss contribution yields 

fp = \/2/p,em fm = 3y/2f wem f§ = —3fij),em 

where (Q.\J^\V,p,e) = myfyE^p) defines the isospin projected decay constants and the current is 
understood to correspond to the given meson. The resulting numerical values are 

f p =221.3(1.0) MeV /„ = 194.7 (3.2) MeV / = 228.7 (3.6) MeV , 

which are all of the same order of magnitude. This is the expectation from chiral perturbation 
theory assuming precisely this pattern of (0-(j) mixing. Letting fy denote the generic vector meson 
decay constant in this case gives 

f — Jx- f — f v f — -ix. 

This can be generalized to include the ft)-0 mixing with an angle 6, resulting in 
f - f y f ■ a f y f a f y 

fp,em — —^ fm.em — sm "~^ fij),em — COS K)—j=^ . 

As needed shortly, we note that tan 6 = 1/ y/2 is the angle that corresponds to a pure ss state for the 
and this angle reduces the previous line to the line before. 

3.2 Vector meson contributions to Aa M 

The couplings in the previous section are on-shell properties of the mesons, but now we must 
discuss the off-shell aspects of the vectors. This inevitably introduces a model dependence. Tree- 
level calculations in chiral perturbation theory provide a definite off-shell form of the vector meson 
propagator and give a contribution to the renormalized vacuum polarization function of 



2\ fv,em Q 2 



rriy Q 2 + niy 
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This expression can also be achieved by simply assuming an off-shell propagator of the form 
A/(Q 2 + my) and fixing A by demanding the correct T(V — > e + e~) result. 

Combining the above result with the relationships between the decay constants gives 

n r^-fJ 1 Ql,m2 p i sin2e Ql K i cos20 g2/m ^ 1 
n^(fi)-/, ^+ 6 Q2+< + 6 Q2+m 2 >- 



Additionally setting m p = m m = m§ = my gives 

II v ,<J- - M 



2$ Q 2 
3m 2 Q 2 + m 2 v 



In the above we now make explicit that this is envisioned as a three-flavor result. We can also 
extract the two-flavor result by decoupling the and demanding a pure ss state for the by setting 
tan 6 = 1/ \pl. The resulting expression is 

rr (n 2 ^ — 5 $ 



9m v Q 2 + m 2 v ' 

Notice that the strength of the Nf = 2 and 3 results follows the sum of the charges squared. 

The resulting integral in Eq. 2A_ can be performed giving the following contributions for each 
of the vector mesons. 

Aa M , p = 470.9. (4.0) • 10~ 10 Aa^ a = 39.1 (1.3) • 10~ 10 Aa M , = 39.0 (1.3) • 10~ 10 . 

These values already reproduce much of the experimentally determined contributions given earlier. 
The discrepancy is largest for the p and is a likely indicator of the effects of the p — > Tin decay. 



4. Lattice calculation 

The vector decay constant can be calculated from the correlator of the non-singlet current as 
follows. 

£ jd 3 x(J ! =\t,x)jl=\0)) -> t^^ {m vMV,e\m) = l^e- m * 

This correlator can be calculated on the lattice without disconnected diagrams. Let C jlv (x) de- 
note the connected piece of the single quark correlator (7^(^:)7v(0)). Similarly let D llv {x) be the 
disconnected piece. Then the 1=1 correlator is 

(J l =\x)j'= 1 (0))=C, v (x). 

Additionally, the electromagnetic coupling of the p can also be calculated directly without discon- 
nected diagrams. 

(jfr(x)J i = i (o)) = ^=c, v (x) 

This is consistent with the earlier relationship f Piem = f p /yj2. 
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L/a 


a}X 


m n [MeV] 


m v [MeV] 


fv [MeV] 


Aa^y [10- 10 ] 


3.90 


20 


0.0040 


347.9 (6.2) 


1209. (161.) 


345. (64.) 


209. (45.) 


3.90 


24 


0.0150 


645.9 (1.7) 


1235. (23.) 


322. (10.) 


164. (6.) 


3.90 


24 


0.0100 


524.5 (1.2) 


1156. (36.) 


312. (13.) 


199. (11.) 


3.90 


24 


0.0085 


484.6 (1.2) 


1144. (34.) 


305. (13.) 


198. (9.) 


3.90 


24 


0.0064 


423.1 (1.0) 


1083. (35.) 


297. (11.) 


234. (15.) 


3.90 


24 


0.0040 


340.2 (1.7) 


1067. (46.) 


290. (15.) 


237. (19.) 


3.90 


32 


0.0040 


334.2 (0.5) 


1044. (51.) 


296. (17.) 


267. (26.) 


3.90 


32 


0.0030 


291.5 (1.0) 


956. (70.) 


260. (20.) 


295. (45.) 


4.05 


24 


0.0060 


453.5 (3.4) 


1094. (39.) 


308. (14.) 


241. (15.) 


4.05 


32 


0.0080 


517.1 (1.6) 


1124. (36.) 


298. (12.) 


204. (13.) 


4.05 


32 


0.0060 


448.5 (1.9) 


1104. (45.) 


299. (17.) 


219. (14.) 


4.05 


32 


0.0030 


325.1 (1.9) 


1063. (78.) 


297. (25.) 


254. (35.) 



Table 1: Ensembles and results used in this work. The pion masses are taken from [10] with one exception. 
The L/a = 20 result was communicated privately by the ETMC collaboration. The lattice spacings are 
a = 0.079 (2) fm and 0.063 (1) fm for J3 = 3.90 and 4.05 [10]. Results from this work are the vector meson 
masses ray, decay constants fy and the corresponding contribution to Aa^ (the last two are normalized as the 
p contribution). Up to disconnected contributions, the Nf = 2 value of Aa^ is 10/9 of Aa^y. Additionally 
ignoring s quark quenching effects, the Nf = 3 value is 4/3 of Aa^y. 

The ft) and <p decay constants are then obtained from 

(^=°(x)4=°(0)) = C MV (x) + 2D MV (x) 

and 

(^(x)4(0))=C MV (x)+D MV (x). 

In the last expression we are assuming a quenched strange quark that is degenerate with the u and 
d quarks. Thus ignoring disconnected diagrams we find m m = m§ = m p and f m = f§ = f p . The 
electromagnetic couplings can also be calculated using the following expressions. 

</-(*)4=°(o)> = ^ c„ v (*) + 

(7-(x)4(0)) = -ic MV (x) 

Again we see that up to s quark quenching effects and disconnected diagrams the results f a ,em = 
fco/ (3\/2) and f§, em = — /3 are reproduced. This makes it clear that the theory with a quenched 
strange quark should be compared to the tan 6 = 1/ \pl scenario in which the is pure ss state. 
The full electromagnetic vacuum polarization function in the Nf = 2 calculation is 

<7-(x)/n0)> \ Nf=2 = \c, v {x) + l -D, v {x) 
and the result including a quenched and degenerate strange quark is 

(jfr( X )J7(o))\ Nf=3 = lc, v (x). 
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Figure 1: Vector coupling. The vector coupling, nor- 
malized as the p, calculated for a = 0.079 fm and 
a = 0.063 fm agree. A linear extrapolation in irv\ 
to the physical pion mass agrees with the value ex- 
tracted from r(V — > e + e~). The factor of my in fy 
hides this agreement. 



Figure 2: Vector contribution to Aa^. A quadratic 
extrapolation in m\ undershoots the measured contri- 
bution of the p significantly. A model (described in 
text) suggests the importance of calculating at small 
enough m K for the vector mass to approach the phys- 
ical p mass. 



Again notice that, up to disconnected contributions and quenching effects, the Nf = 2 and 3 results 
agree to within the sum of charges squared. 

5. Results 

Our earlier calculation [6] has been extended to include a second lattice spacing, more quark 
masses, volumes and estimates of disconnected diagrams. (See Tab. |].) Here we focus on the sim- 
plest model for Hr(Q 2 ) that includes only the vector mesons. The vector meson masses and decay 
constants and the corresponding contribution to Aa^ are given in Tab. [I[ The vector contribution 
to Aa^, normalized as the p, rises from 164 (6) • 10 10 to 295 (45) • 10 10 but falls well short of 
the expected 501.1 (3.0) • 10~ 10 . 

We suspect that this discrepancy is mainly due to the large value of the p mass found for the 
range of n masses in Tab. [l[ To illustrate this point, we plot the coupling of the vector meson 
in Fig. U There we see that gy appears to have a mild quark mass dependence and a simple 
linear extrapolation in m\ gives a value in agreement with the experimental measurement. We can 
isolate the dependence of Aa^y on gy and my by integrating Eq. ^j] exactly and then expanding 
in m^i/my. The first few terms are 

\my J I 3 12 \my J \my J \my J \my J J 

This makes it plausible that the discrepancy is simply due to the large value of the p mass and 
the approach to the physical point should account for most of the discrepancy between the current 
values of Aa^y and the experimental measurement. To illustrate this, we fit the chiral expansion 
my = a + bm 2 l +cm\ to both our lattice results and also the physical value of the p mass. Combining 
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this with the fit to gy we produce the model curve in Fig. ^[ Because gy extrapolates to the physical 
value and my was constrained to give the physical p mass, this model automatically reproduces the 
narrow p approximation from earlier. 

6. Conclusions 

The final results of this work will be given in a forthcoming publication. Here we focus on 
understanding the dominant contribution from the vector mesons. Our current calculations produce 
values of Aa^ that are significantly lower than the experimental measurements. We construct a 
model that suggests that this effect may be caused partially by a p mass that is still rather large 
compared to its physical value. Thus we expect that calculations with the physical value of the 
71 mass and including the strange and charm quark contributions will be capable of achieving the 
precision necessary to match the accuracy of current experimental measurements of Aa^. 
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